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Local Heavy Quasiparticle in Four-Level Kondo Model 
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An impurity four-level Kondo model, in which an ion is tunneling among 4-stable points 
and interacting with surrounding conduction electrons, is investigated using both perturbative 
and numerical renormalization group methods. The results of numerical renormalization group 
studies show that it is possible to construct the ground state wavefunction including the excited 
ion states if we take into account the interaction between the conduction electrons and the 
ion. The resultant effective mass of quasiparticles is moderately enhanced. This result offers a 
good explanation for the enhanced and magnetically robust Sommerfeld coefficient observed in 
SmOs4Sbi2, some other fiUed-skutterudites, and clathrate compounds. 
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1. Introduction 

Since the two-level Kondo effect was studied about a 
three decades ago,^'^ the physics arising from the off- 
center degrees of freedom of ions have attracted much at- 
tention. Recently, off-center motions of ions in clathrate 
compounds'^ and filled-skutterudites^ are identified by 
ultrasonics. A dispersion of the sound velocity and a soft- 
ening in the elastic constants are observed in various ma- 
terials such as PrOs4Sbi2, LaOs4Sbi2, Ce3Pd2oGe6, and 
so on. Characteristic properties of these materials are 
the existence of a thermally excited rattling motion in 
the "cage" , and a possible realization of an anisotropic 
ground state with respect to the ion configurations. For 
example, the elastic constant C44 of La3Pd2oGe6'* dis- 
plays a softening at least down to 20mK. This implies 
the existence of a triply degenerate mode in the ground 
state not in the singlet state which is expected in a simple 
quantum mechanical system. 

About two decades ago, the energy level scheme for 
a simple and canonical off-center mode was obtained in 
NaCl including 0H~ in its interstitials.^' ^ The OH^ ion 
is located at the [100] off-center positions in NaCl. The 
result of the absorption spectrum of a millimeter wave 
supports a singlet ground state (T^) with an excited 
state triplet ). In that case, no softenings in the elas- 
tic constants is observed at low temperatures, reflecting 
the singlet ground state. 

The off-center positions in PrOs4Sbi2 and LaOs4Sbi2 
are expected to be the same. In the case of PrOs4Sbi2, 
in addition to the f-electron contributions, the elastic 
constant Cn — C12 continues to display an anomalous 
softening to the temperature just above the supercon- 
ducting transition temperature.^ This implies that, in 
a low temperature region, there exists a r3(r23) mode 
of the ion configurations, while an elementary quantum 
mechanics tells us that the lowest energy state is the 
singlet with a symmetric ion configuration in such tun- 
neling systems. In this point, the observed results are 
mysterious. It is natural to expect that a magnetically 
robust heavy fermion compound SmOs4Sbi2^ is added 



into such a group of materials. 

As theoretical developments for the two-level systems 
(TLS), it was pointed out that a strong interaction be- 
tween local phonons and conduction electrons gener- 
ates such off-center situations as a model of A15 com- 
pounds. ^"^^ Recently these problems are also discussed 
in the lattice model in relation to the rattling degrees of 
freedom in fiUed-skutterudite compounds. The prob- 
lems of the off-center motions of ions in the two-level 
system coupled with conduction electrons were studied 
extensively^^"^^ in the context of the two-channel Kondo 
model. 

The purpose of this paper is to study a property of 
a four-level Kondo model which is a simplified model 
for the fiUed-skutterudite compounds with apparent off- 
center degrees of freedom of ions in the cage. In particu- 
lar, a condition for an appearance of heavy quasiparticles 
is investigated on the basis of not only the perturbative 
renormalization group (PRG) method but also the nu- 
merical renormalization group (NRG) one. 

The systems with off-center positions more than two 
were also discussed so far. A/-level system was discussed 
in ref. 20 and a symmetric SU(M)(g)SU(A^/) Coqblin- 
Schrieffer type effective Hamiltonian was derived using 
expansion, where Nf is a number of the conduction 
electron channels (spin). In Ref. 20, the only conduction 
electron orbital which couples most strongly to the ion 
is retained. In this paper we investigate the properties, 
which arise from the purely orbital degrees of freedom of 
the conduction electrons. 

This paper is organized as follows: In § 2, we derive 
the effective model Hamiltonian for a four-level Kondo 
model. The details of explicit derivation of four-level 
Kondo model are given in the Appendix A. In § 3, we 
discuss the property of the four-level Kondo model using 
scaling approach on the basis of PRG. In § 4, we discuss 
the four-level Kondo model using NRG method. In § 5, 
we give a possible explanation for the unusual properties 
observed in the above-mentioned materials. Finally, § 6 
summarizes the result of the paper. 
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Fig. 1. (a, 0, 0) type off-center positions of an ion. Each circle 
represents off-center positions located at (±a, 0, 0), (0, ±a, 0) and 
(0,0, ±a). The numbers correspond to indices of 0i(x). (a)two- 
dimensional square case, (b) three-dimensional cubic case. 



its Fourier transform and the system volume, respec- 
tively. F^j (q) is given by 



d3:re-"i-(/)f (x)0°(x) 



(6) 



■1 exp[-^ — 1^ — ^ • q] 



4a2 



(7) 



Before writing a Haniiltonian in this system, we should 
orthogonalize ^^(x).^^ To leading order in S, we obtain 
the orthogonalized wavefunctions: 



.(x)^0°(x)-i5^0°(x). 



(8) 



2. Models 

In this section, we derive an effective model Hamil- 
tonian describing the four-level off-center mode coupled 
with conduction electrons. 

2.1 Effective model Hamiltonian 

First, we assume that the ground-state wavefunction 
of an ion localized at each off-center position x^ {i = 
1, 2, 3, • • • ) is approximated by the Gaussian form as, 

\2. 



1 



exp 



(x- 



2a2 



(1) 



where x is the coordinate of the ion, cr is a broaden- 
ing length of ion wavefunction. In this paper, we investi- 
gate the case where the off-center positions are (±a,0), 
(0, ±a), see Fig. 1(a). The state \i) corresponding to 0°(x) 
is mixed with |j) (i ^ j) by the quantum mechanical 
tunneling. Note that we consider the low temperature 
physics so that thermally excited classical tunneling pro- 
cess are discarded, or such processes are renormalized in 
the parameters of the low temperature effective Hamil- 
tonian. The overlap Sij between states |i) and |j) are 
estimated by an usual integral as follows: 



= / d^x 0r(x)0°(x) 



(2) 



52 



= S* <C 1 for the nearest neighbors^g-j 
for the next— nearest neighbors. 



Because of property (3) we retain only the nearest- 
neighbor hopping for simplicity. These tunneling terms 
cause the energy splitting with respect to the irreducible 
representations of ion wavefunctions, as discussed in the 
next section. 

Next, we consider conduction-electron assisted tunnel- 
ing. These processes are occurred by the following matrix 
elements, 



d^a;d^r0r(x)e~"'""t/(x - r)e"''""0°(x) 



(4) 



= ^;°(k-k');7k-k'. (5) 

Here, C/(x — r), ?7k-k' and are the screened Coulomb 
interaction between the ion and the conduction electrons. 



Here, the summation of j runs over the nearest-neighbor 
off-center positions of i. 

From eqs. (3) and (7) the Hamiltonian describing the 
off-center motion of the ion and the interaction between 
the ion and the conduction electrons are given as 



H 



fkCk^Ck^ + Hint , (9) 



H^nt = 5]C/k-k'4ck>5]F,,(k-k')40.(lO) 

kk'/i ij 

where c^^ and 0i is the annihilation operator of the con- 
duction electron with the momentum k, spin /i and that 
of the ion (orthogonalized) with the off-center position 
Xi, respectively. The tunneling amplitude Ay is esti- 
mated as Ajj ~ —hujoSij where luq is the frequency of the 
ion around each off-center positions. Fij(k — k') is esti- 
mated by using instead of 0°(x) in eq. (6). Note 
that an ion operator (f)i acts on a restricted Fock space 



with J2i 



At. 



1. Non-commutative properties of this 



model come from the off-diagonal terms of -ffint, which 
give rise to, in simple cases, various Kondo effects. 

In the following, we ignore the spin of the conduction 
electrons and expand Ck by the spherical harmonics at 
the origin as in the usual way of an impurity model. Phys- 
ical properties generated from the spin dependence will 
not be important because the resulting effective Hamil- 
tonian is highly anisotropic in the orbital space of the 
conduction electrons (s-, p-, d-wave), e.g., giving little 
chance for the two-channel Kondo effect to be realized. 
Then we expect that criticalities such as the two-channel 
Kondo model cannot arise in the present model. More- 
over, The existence of the splittings Ay also helps this 
expectation. 

2.2 Order estimation of parameters 

Before proceeding to detailed analyses, we need to es- 
timate the order of magnitude of parameters in our effec- 
tive Hamiltonian eq. (9). Using a one-dimensional double 
well potential. 



Vix) 



(11) 



the stable points ±|a;*| are estimated as |a;*| — 1/273 

aw~^, where w is a parameter of 0{1). We can estimate 
the frequency loq of the harmonic oscillator centered at 
±|a;*| as hujo ^ 8wVb, where Vb is the height of the 
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barrier at x ^ 0. For the nearest-neighbor pair 
Ay ~ huoaS = Sw^VbS and [/(q)Fy(q) - [7(q)5 ~ 
AjjC/(q)/(?iwo) = Ayi7(q)/(8VBW^). These expressions 
coincide with the result of the detailed analysis in an 
order of magnitude given in ref. 21 . In the case of 
PrOs4Sbi2, La3Pd2oGe6 and Ce3Pd2oGe6, Vb is esti- 
mated as 70 — 200K.^ Taking into account the fact, 
a ~ |a;*|, 5 ~ exp(-l/(2w;2)) = 0.0183 for w = l/{2^/2). 
Then A^j ~1K in the case that is in the near- 

est neighbor with each other. An ambiguity about the 
interaction term t/(q) remains. There are two limiting 
cases: \U{q)\ > Vb i\U{q)\ ~ lOOOK) and |;7(q)| < 
Vb (|f/(q)| ci: lOK). Detailed investigations will be given 
in later sections. 

3. Perturbative RenormEilization Group Ap- 
proach 

In this section, we calculate the effective interactions 
using PRG method, for the purpose of seeing the overall 
features of the scaling in our model. Although there ex- 
ist some studies by using renormalization group in such 
models, ^'^ our approach is a natural extension of rcfs. 13 
and 14. A change from refs. 13 and 14 is a treatment for 
bandwidth 2Do of the bare conduction electrons. It was 
pointed out in ref. 23 that Dq should be replaced by an 
energy level of the first excited state (es in their notation) 
which is smaller than Dq {e^ ^ (10~^ ~ 10~^) x Dq). 
Because of this, it is concluded that the Kondo tem- 
perature is always smaller than the spontaneous split- 
ting of the ion configuration, in weak coupling. Recently 
a mechanism of a resonant scattering of the tunneling 
impurity was proposed in order to obtain rather high 
Kondo temperature.^^ In this paper, we consider the case 
Do ^ 100K~ Vb- Note that Vb is comparable to Huq 
which is an energy spacing of the local phonon, i.e. Ein- 
stein frequency. 

3.1 Effective exchange Hamiltonian 

We expand the annihilation operator Ck of the conduc- 
tion electron in terms of the spherical harmonic function 

I 

Ck = ^ ^ Yim{k)i''Cklm, (1) 
I m=—l 

where we define a new annihilation operator Ckim with 
k = |k|, the angular momentum / and its z-componcnt 
TO. In eq. (1) we have omitted a coefficient depending on 
k,^^ which causes no significant problem by absorbing it 
in a definition of the bandwidth 2Dq (centered at the 
Fermi energy). The interaction f/k-k' is also expressed 
by the Legendre expansion as, 

f/k-k' =E E (2) 

I 711 = — I 
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where ui is the /-th coofBciont of the Legendre expansion 
for f/k-k'- Then Hamiltonian (9) can be written as 

H = ^A^ata^ + ^efcc;f-Cfc;+ifi„t, (3) 

A* hi 

= Y.Y.Y.T.%ifu'P<n^-- (4) 

fefe' III nv 7 

Here, we have introduced new ion annihilation operators 
a^'s which arc the irreducible representation of the point 
group with the energy level A^ of the ion configuration 
(§ 3.2). The indices I and V represent the angular mo- 
mentum {l,m). The coupling constant J7~ is calculated 
in Appendix A for the four-level system, 'if'^^ is matri- 
ces of the 7-th irreducible representation of the direct 
product aj^a^, and normalized as X^/^y 1^/11;/ P = 1 (§ 3.2). 

For this kind of rather general exchange Hamiltonian 
(4), a formula of 2-loop renormalization group equations 
has been derived in ref. 25. The renormalization group 
equation for A^ is obtained in a way similar to the case 
of TLS. Then a set of renormalization group equations 
is given as follows: 

-^EE-^^'4-'T<[*".*^]*^0 

l" a/3 

+ EyTr(j"J'')jATr([*",[*^',*^]]*^t^, (5) 

^= vet<^"^'^)e*;^^*^.^- (6) 

a/3 ^^T^A* 

where x = log(I?/Z?o), D being the half of the scaled 
conduction-electron bandwidth, and p is the density of 
states (DOS) at the Fermi level. We take the same band- 
width and the DOS for the all partial-wave components 
of the conduction electrons. The renormalization group 
equations, (5) and (6), are quite general and are not re- 
stricted in details of models. In the following, we apply 
these general expressions to the four-level Kondo model. 

3.2 2- dimensional m,odel 

In 2-dimensional configuration of Fig. 1(a), we trans- 
form the ion operators 0,'s to the irreducible representa- 
tions a^jjs of Dih point group as 

r+: Oo = \[(l>i+h + (t>3 + (t>il (7) 
r+: 02 = \['^i-h + <t>z-(t)i], (8) 

r^± : a±i = ^[(^1 T#2 - </>3 ±#4]- (9) 

Furthermore, we need to define the operator form of "^"^ 
(which is defined as = aj^'^^^o,^) for the direct prod- 
uct aj^tti,: 

2*° = alao + a\ai + ala2 + a^_-^^a-i, (10) 
\/2*^ = ajao - 4a2, (11) 
\/2*^ = a\ai - aLifl-i, (12) 
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2*^ 




a^ao — a[ai + w^ai — a_]^a_i, 


(13) 


2*^ 




J. J. J. J. 
ajao + a^ai + al_ja2 + Oq'^-i, 


(14) 


2*^ 




iff f 

a|ao — 02*^1 — a_^a2 + a\^a-\^ 


(15) 


2*^ 


= 


— ajoo + a2'^4 — 0,1^^0,2 + a\a^i. 


(16) 


2*^ 


— 


— ajao — 0301 + al]^a2 + agO-i, 


(17) 


^8 


= 




(18) 




= 




(19) 




= 


[*«]t, 


(20) 








(21) 


\/2*i2 




a|a_i + al^j^ai, 


(22) 


\/2*i3 




a|a_i — a]^]^ai, 


(23) 


x/2*i^ 




aja2 — 012^05 


(24) 


\/2*i^ 




aja2 + 02^05 


(25) 



with constraint 2#° = 1. For the conduction electrons, 

we truncate the summation with respect to Z at Z = 2 in 
eq. (1), because dominant contributions are expected to 
arise from the electrons with low I waves. Corresponding 
irreducible representations for the quadratic form made 
from these conduction electron field operators are listed 
in Table I. 



rep. (-D4h) 


conduction electron 


ion 




sts, 

stdg + h.c, 
d^^d+ + did-, 
+ h.c, 

P+P+ + vlv- 






p\p+ - pip-, 
— dl (i_ , 
(i^_d_ — h.c. 


^2 


^ 3 


+ h.c, 
s\d^ + c(_) + h.c, 
(ij(d+ + d-) + h.c 


^12^15 


^ 3 


st(d+ + d-) - h.c, 

dj{d+ + d_) - h.c. 




n 


p^p_ — h.c, 
st(d^ — d_) — h.c, 
dQ(d+ — d_) — h.c. 


*13 


^ 5 


{-stp+ + pis, stp_ - pis}, 
{-djp+ +pldo,dJp_ -p^i^o}, 
{c^lp- -pld-,-dlp+ +pld+}, 
{dlp_ — pYd+, -djj_p+ +pld_} 


1*8^^4} 


r(2)- 
^ 5 


{— stp^ — pls,stp— +p^s}, 

{dtp_ +pU-,dlp++pld+}, 
{dip- +p+d+,d^,_p+ +pld_} 





Tabic 1. Irreducible representations for the quadratic terms of 
the conduction and the ion operators. We classify them by tak- 
ing into account the fact that the couplings are all real quan- 
tities. In the second column, we use abbreviations such as s = 
Efe Cfeoo, P± = Efc Cfciii, do = Yjk ^fc20 and d± = Yl,k c/c2±2- 



Then, the "bare" Hamiltonian is (see Appendix A for 
detailed derivation): 



^diag H~ ^hopi 
2kFCL 



(2 



k'i.a^ , t t 
^{p+P- +p'-P+ 



-{s^{d+ + d-)+ h.c. 



*12 + * 



15 



x/2 



(2 



op 



^^M^{2st5+^(stdo + 4s) 
o v5 



\/2k%a'^ f I + 

- P-P+ 



d_)- h.c. )}*i3, (28) 



where we have used the abbreviations in Table I for the 
conduction electron field operators, and keep only I = 
term in eq. (2). In cq. (26), there exist no Ckia and Ck2±i 
terms which emerge in the higher order interaction, e.g. 
in the case when we consider wi . However, these terms are 
expected to be less relevant because of the 2-dimensional 
configurations of the ion (i.e., the ion is located on the 
xy- plane). For the local ion energy, we obtain: 



A(ojao - 4a2) (A < 0). (29) 



In Fig. 2, we show the result of the 2-loop rcnormal- 
ization group flow for a typical case. It is noted that the 
1-loop result is different from 2-loop one especially in 
•^p+d+ ^'^'1 Jd+d+- The fixed point of 2-loop calculation 
is an artifact arising from the third order approximation 
in J2, as in the single-channel Kondo effect. Detailed 
studies show that the fixed point given by a 1-loop 
result might be realized. At the fixed point of the 1-loop 
calculation (and of course around the same energy re- 
gion for the 2-loop one), dominant coupling constants 
are those between (p_)_,s,p_) and (ai,ao,a_i), and be- 
come strong coupling. 

For 7io > 0, the effective Hamiltonian can be written 
in the following form: 



1 



iir« = Ji(s,5, + ^(gf oL') + ^)gf))-hJ2(^g.Q. 

S. + S+S+ + ^I'^OL'^ + q^'^Qf)), (30) 



where J2 — jJi > 0, and Sz{Sz) and s±{S±) are the 
z- and transverse-component of pseudospin s{S) = 1 for 
the conduction electron (ion). The eigen states of Sz and 
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Fig. 2. 1- and 2-loop renormalization flow of the coupUng J^^ 
(only typical couplings are shown). For example, a coupling 
•^pXp— niGans one corresponding to the term p^p_^^^. The 
Kondo temperature — 0.009-Do of 1-loop result is shown 
by a vertical dotted line. The parameter set used is: uq = 
0.6Do,kpa = 0.8, kpa = 0.42, and A = -O.OIDq. The val- 
ues of J^j's below ~ Tk or have no meaning because PRG 
would not be valid there. 



Uo=M 




D/Dr 




10 T/D„ 10 



Fig. 4. Cimp/T vs. T, C being the impurity specific heat. The 
parameter set used is: A = — 0.0169So, kpa = 0.8 and 
kpa = 0.42 (we obtain similar results for jiq > cases). 
For the bare bandwidth Do/ks = lOOK, d^p/T ~ 57 X 
(scale of ordinate) mJ/K^ ■ mol. 
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D 
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a 
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Uo/Do= 



•-•-•-•-•-•-•-•-g;0=o-o-o;2;a;a=fl 

O-O-O-O-O-O-O'O' kA 



10"' 10"' 10"' io't/DqIo" 



Fig. 5. Occupation number of ionic configurations F's as a func- 
tion of temperature. The parameters used are the same as those 
in Fig. 4. Although we show results for mq < cases, we obtain 
similar results for «o > cases. 



Fig. 3. Renormalization evolution of the energy splitting of the 
ion configurations A^ as a function of the scaled bandwidth 
D/Dq. Arrows indicate the Kondo temperature Tk given by the 
1-loop calculation. The calculated A^^'s have no meaning below 
~ Tk or A^. 



Sz are constructed as 

|s ===1,5,-0) = st|o), (31) 

|s = l,s, = ±1) = ±p^±|0), (32) 

1^=1,^, =0) = a>), (33) 

\S=l,Sz = ±l) = alilO). (34) 

On this basis set, s and q are expressed by 3 x 3 matrices 
as follows 

1 \ /O V2 

,s+: \/2|, (35) 

00-1/ \0 










2 








73 , 



V2 


















,f : I 



and S- 





(1)^' J (2) 

and q_ 



(2)t 



(37) 



The operators 



g's represent the quadrupolar degrees of freedom. Sz,Qz, 
etc. are defined in the same way. In these representations 
the ion energy of eq. (29) implies a local distortion field 
in pseudospin space, because the following relation holds: 



ala2 



VSQz 



(38) 
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1 2 3VDo 

Fig. 6. Phase shifts of the conduction electrons with t;- 
symmetry vs. the interaction uq. The parameters used are the 
same as those in Fig. 4. A similar result is obtained for ug < 0. 
The small splitting between and p_ at uq = SDq is due to a 
numerical error from truncation procedures in NRG. 

near the strong coupling fixed point where 0302 is irrel- 
evant (the degrees of freedom are frozen out in low 
temperature). 

For Mo < Oj we also obtain an expression for the effec- 
tive Hamiltonian similar to eq. (30). From eq. (30), it is 
expected that an anisotropic spin 1 dipolar {sz , s± ) and 
quadrupolar (g^, q± \ q±'^) Kondo effect with the conduc- 
tion electron with the pseudospin 1 occurs in this model. 
For comparison, the more symmetric Kondo model, in 
which the impurity and the conduction electron with 
spin 1 are coupled isotropically, has two fixed points. 
One is the SU(3) symmetric fixed point, which is oc- 
curred when we take into account both the dipolar and 
the quadrupolar interactions. The other is a non-Fermi 
liquid one keeping the only dipolar interaction. 

In Fig. 3, the energy levels of the various ion config- 
urations are shown. A renormalization of does not 
occur at the level of 1-loop PRG (poorman's scaling) 
but appears at the 2-loop PRG equation eq. (6). It is 
noted that the energy level A±i of the first excited state 
decreases as the effective conduction electron band- 
width decreases. On the other hand, that of the sec- 
ond excited state Fg , A2, increases. The results in the 
low D/Dq region are unphysical ones, because the per- 
turbative approach cannot be reliable there (especially 
i5<max(|A^-Ao|,Tx)). 

4. NRG Results 

In this section, we investigate the 2-dimensional model 
(eq. (26)) using the NRG method. Because of compu- 
tational restrictions, we take the only s, p± and d± con- 
duction electrons into account. The omitted orbital do is 
expected to play an irrelevant role from the result due to 
the PRG study discussed in the previous section. For the 
NRG calculations, it is helpful to introduce the following 
two conserved quantities: 

Q ^ E(4i™CM™ - 1/2), (1) 

klm 
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h = modC^mcli^Ckim+ M^'^M' 4), (2) 

klm M=±l,0,2 

where Q is the total charge of the conduction elec- 
trons measured from the half filling and the helicity h 
is a kind of the z-component of the angular momentum. 
In the NRG calculation, we transform the conduction- 
electron operators Ckim into the 1-dimensional semi- 
infinite chains fnim {n — 0,1,2,...) and diagonalize it- 
eratively the Hamiltonian starting from the n = shell. 
The conduction-electron part in the Hamiltonian (9) is 
expressed as follows: 
00 

^kcli^Cklm = tn^^"^ fli,Jn+llm + h.C, (3) 

klm n—0 Im 

t = (i + A'^)(i-A— ^) 

2^(1- A-2«-i)(l- A-2«-3)' 

where A is a logarithmic discretization parameter in the 
NRG calculation. We keep 1200 states at each iteration 
and use Dq = 2Dq/{ \ + A^-'^) as a unit of the energy and 
set A = 3 in this section. 

In Fig. 4, we show the result of Cimp/T the impu- 
rity specific heat (divided by temperature). Although we 
show results for mq < case in Fig. 4, we obtain simi- 
lar results for uo > case. In the region of small mq, we 
do not have a heavy effective mass, because the Tk in 
this region cannot exceed the local splitting A. In the 
region of large uq, say for uq/Dq = —3 in Fig. 4, the 
Tk can overwhelm the splitting A. In this case, however, 
the Tk is too large so that the effective mass becomes 
small. As a result, we have the heaviest effective mass 
at an intermediate strength of uq i.e., uq/Dq = —2 in 
Fig. 4. These features are common in the usual Kondo 
model under a magnetic field, which has an equivalent 
function to the local splitting in the present model. A 
size of Ci„ip/T in a physical unit depends on how to 
take Do{Do). If we set Do/kB — lOOK as discussed at 
the beginning of § 3, it is given as Cimp/T ~ 57 x 
(scale of ordinate) mJ/(K'^-mol). Then for uq — —2Do, 
limr^o Cnnp/T ~ 4.6 x 10^ mJ/(K2-mol), which implies 
that local heavy electron state is realized. 

The occupation number for each configuration of the 
ion is shown in Fig. 5. As expected from the results in 
the last section, a Kondo effect between (oi, oq, a_i) and 
{p+, s,p-) seems to be occurred. It is clearly seen in Fig. 
5 that the orbitals, which take part in the Kondo ef- 
fect, grows at low temperature. It is noted that the state 
r^j_, originally doublet excited states, are included in the 
manybody ground state together with F^'" fully symmet- 
ric state. This is marked contrast with a simple picture 
of quantum tunneling among multistable points and also 
with the case of two-level Kondo problem in which an 
ionic configuration of the ground state is fully symmet- 
ric like a simple TLS. 

In Fig. 6, we show the phase shift of conduction elec- 
trons (^ = s,p±,d±). We estimate using the energy 
level of the first excited states in each ^. It is noted that 
Sd± ~ for any values of uq while (5^ (^ = s,p±) are 
enhanced as uq is increased. This confirms the result by 
the PRG approach discussed in § 3, where it has been 
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shown that d± orbitals do not take part in the fixed point 
Hamiltonian and are decoupled. 

5. Discussions 

As discussed above, a realization of Kondo screening 
in multilevel systems crucially depends on the strength 
of the coupling constant between the conduction elec- 
trons and the ion charge fluctuations. In the last section, 
we have used the couphng mq = 0(1) x Do- It appears 
that such a coupling is unrealistically large. However, we 
have to remember that Dq is an effective bandwidth of 
the order hwo. As argued by Aleiner et al.^^ the conduc- 
tion electrons with energy ^ fvjjQ may not contribute a 
renormalization of the coupling ui in the present model 
when the effective bandwidth is scaled down. As a re- 
sult, we can use relatively strong coupling constants and 
splittings of ion's level. Since the spin degrees of freedom 
have been neglected in the present model, the enhanced 
effective mass obtained here is evidently robiist against 
a magnetic field. It is interesting to note that the heavy 
electron state, which is robust against the magnetic field 
up to 14T, has been observed in SmOs4Sbi2.^ 

Although the present model (2-dimensional version) 
cannot be compared directly to the real materials, 
the important result is that it is possible to include 
anisotropic ion modes in the ground state wavefunction 
if the coupling between conduction electrons and the ion 
charge fluctuation is moderately strong. This "inclusion" 
is expected even at the level of 1-loop renormalization 
group. Indeed, since the effective Hamiltonian at the 
fixed point is written as an anisotropic spin-1 Kondo 
Hamiltonian, the selected ion's modes and Fg^. have 
energy gains by the Kondo effect. In 3-dimensional model 
(Fig. 1(b)), a similar renormalization is expected to occur, 
involving the anisotropic F4 triplet configuration (excited 
state) together with the isotropic Fi configuration. 

In order to discuss the thermal rattling, it becomes 
necessary to include more realistic aspects in the model, 
e.g., the role of excited states of the ionic motion and 
the effect of dissipation in the local system. This is be- 
cause, in ultrasonic experiments, the frequency u> of the 
ultrasonic is a; ~ ImK, so that it is difficult to realize 
the dispersion of clastic constants at around 30K with 
respect to the frequency"^"^ only by taking into account 
the complete energy levels of the local system without 
dissipation. An important ingredient giving the disper- 
sion would be the dissipation due to the interaction be- 
tween the electron and the local ion, or that between the 
local ion and the acoustic phonons. Indeed, in refs. 3-5, 
a Debye-type dispersion is used to analyze the results of 
the experiments. The reason why the ultrasonic disper- 
sion (rattling) is observed only in specific modes might 
be related to the matrix elements of the dissipative in- 
teraction. 

It is noted that there are two aspects in this prob- 
lem. One is the physics of the heavy fermion in which 
the several lowest-energy levels of the local system and 
the orbital Kondo effect play an important role as dis- 
cussed in this paper. The other is that of the thermal 
excitations coupled to the time-dependent external per- 
turbation such as the ultrasonic, which cause so-called 



rattling. 

For treating excited states, it is needed to solve the 
3-dimensional Schrodinger equation directly, assuming a 
proper off-center potential. Recently, a band calculation 
was performed modeling /3-pyrochlore compounds. The 
eigen value of the energy and the wavefunction for low- 
lying ionic states were calculated, although the couplings 
between the conduction electrons and the ion have not 
been estimated, because the estimation of the coupling 
constants is highly nontrivial. 

In the filled-skutterudite compounds, anomalous 
phonon contributions are observed in various quantities, 
especially in ROs4Sbi2 (R=Pr, La, Sm, etc.). One of the 
reasons why the effect is prominant in ROs4Sbi2 may be 
that the size of Sbi2 cage is the largest among existing 
filled-skutterudite compounds including rare earth ion R. 
Another one may be that the conduction electrons near 
the Fermi level in Os-compounds consist of a molecu- 
lar orbital of the cage (Ai„) and the d-clcctron in Os. 
Unlike Ru-compound, in which the levels of 4d electrons 
are deep compared to that of 5d or 3d electrons so that 
their components near the Fermi level arc; negligible, the 
contributions from d-electrons are much larger. Then it 
is expected that d-electrons play important roles in Os 
compounds. For PrOs4Sbi2. it might be speculated that 
something related to the off-center degrees of freedom 
take an important role for the superconducting transi- 
tion. But the relation between these local phonons and 
unconventional superconductivity is beyond a scope of 
the present paper, and further studies are required. 

6. Conclusions 

In summary, we have studied the four-level system 
with spinless conduction electrons by using renormal- 
ization group methods. The relevant component of the 
ion configurations are Tf and F^ which couple to the s 
and p± conduction electrons in the low temperature fixed 
point. The low energy effective Hamiltonian has been 
written in a pseudospin-1 language both for the ion and 
the conduction electrons with a pseudo-distortion field. 
As a result, anisotropic components of the ion off-center 
modes enter the groimd state wavefunction. We expect 
that the observed relatively large effective mass and the 
robustness to the magnetic field in SmOs4Sbi2 and other 
heavy-fermion-like materials without f-electrons can be 
explained. 

Acknowledgment 

We would like to thank Y. Nemoto and T. Goto 
for fruitful discussions. One of us (K. H.) is supported 
by Research Fellowships of JSPS for Young Scientists. 
This work is supported by a Grant-in- Aid for Scientific 
Research (No. 16340103), 21st Century COE Program 
(G18) by Japan Society for the Promotion of Science, 
and a Grant-in-Aid for Scientific Research in Priority 
Areas (No. 16037209) by MEXT. 

Appendix: A: Derivation of effective Hamilto- 
nian for D4,h symmetry 

In this appendix, we derives the bare Hamiltonian eqs. 
(27) and (28). First of all, Hamiltonian eq. (10) is trans- 
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formed into the form of eq. (4). 

k k' Im I'm' ij 



J J An An 



rm{k)u^F,3{^)yi'nAk') 



•(A-1) 



where q = (k — k') /kp and notice that involves only 
a summation of the radial part k. 

Fij{q) for j m eqs. (10) and (A-1) is approximated 

as 

a? Akpaf 



-Q.Qy + hill + Ql) for = (1,2), (3,4) 



( -Qx 



cQv + Ulx + for (iJ) = (2,3), (4, 1) / '^^'^^ 
where we have expanded -Fjj(q) with respect to kFa{<^ 

2| |2 

1) and approximated exp(— '^"j^' ) ~ 1. In the same 
way, for i = j, we obtain 



1 - i{kFa)qx - 



ql for(i,i) = (l,l) 



1 - \{kFa)qy - ^ql for (i, j) = (2, 2) 

1 + i{kFa)q, - fo^ ^ 3) 



^ 1 + \{kFa)qy 



(kpa) 2 
2 



for = (4,4) 



(A-3) 



In the case of only retaining uq for [/q, the next task is 
to estimate matrix elements oi{qx^qy^ ql^Qxqy, ql) in the 
square brackets of cq. (A-1). For the calculation of these 
matrix elements, it is useful to introduce the following 
spherical harmonics representations: 



X 
^2 



f (n-iTni), 



(A-4) 



3\/5 



Y20±\ —{Y22 + Y2-2) 



-ixy 



§(^2-2 



Y22). 



(A-5) 
(A-6) 



As the matrix elements of 
■ {I'-D f f '^^k d^k' 



"-"// 



An Aw 

we obtain the following results, 

— i 



-Yi*m{k){qx, Qy, ql, qxQy, ql)Yi',m'ik') 



qx 



qy 



AttVQ 

1 

47rV6 



SiiSi>oiSm-i - ^mi) + perm., (A-7) 
i'^;'o(^m-i + Smi) + perm., (A-8) 



+ perm. 



iZn 



-1 T ^ml)(<^m'-l T Sm'l), (A-9) 



qxqy 



-1 



47r\/30 
i 



Si'oSi2{Sm-2 - <5m2) - perm. 



Sll'5li{Sm,lSm'-l - Sm-lSm'l)- (A-10) 

iZn 

Here, "perm." means Im I'm'. Finally, we use eqs. 
(7)- (25), and neglect simple potential scattering terms 
corresponding to the first term in eq. (A-3), then we ob- 
tain eqs. (27) and (28). 
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Note added: 

Recently, Matsuhira et al. [19aYG-13, JPS topical 
meeting, 19-22 September, Tanabe, Japan] reported a 
result of specific heat measurements in LaOs4Pi2- They 
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find that an anomalous increase in the specific heat 
(C/T cx r~2) at low temperatures (T < 2K). Surpris- 
ingly, this behavior is robust against magnetic fields so 
that they claimed it is an intrinsic effect other than the 
nuclear spin contribution. We believe that this power- 
law behavior is explained by the peak structure of the 
specific heat in Kondo systems. Although it is not well 
recognized, Cimp/T in the high temperature side of this 
peak is well described by T^^ law. Indeed, the result of 
Cimp/T by NRG calculation exhibits such T-dependence 
as shown in the figure below. From our point view, the 
anomalous temperature dependence might arise from a 
characteristic of this system in which the very low Kondo 
temperature of the ionic off-center motions is realized due 
to a relatively weak-coupling between multilevel system 
and conduction electrons in P12 cage with smaller sized 
compared to Sbi2 cage. 
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